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ABSTRACT 

Spectroscopic and spectropolarimetric observations of the pre-main sequence early-G 
star HD 141943 were obtained at four observing epochs (in 2006, 2007, 2009 and 2010). 
The observations were undertaken at the 3.9-ni Anglo-Australian Telescope using the 
UCLES echelle spectrograph and the SEMPOL spectropolarimeter visitor instrument. 
Brightness and surface magnetic field topologies were reconstructed for the star using 
the technique of least-squares deconvolution to increase the signal-to-noise of the data. 

The reconstructed brightness maps show that HD 141943 had a weak polar spot 
and a significant amount of low latitude features, with little change in the latitude 
distribution of the spots over the 4 years of observations. The surface magnetic field 
was reconstructed at three of the epochs from a high order (Z < 30) spherical har- 
monic expansion of the spectropolarimetric observations. The reconstructed magnetic 
topologies show that in 2007 and 2010 the surface magnetic field was reasonably bal- 
anced between poloidal and toroidal components. However we find tentative evidence 
of a change in the poloidal/toroidal ratio in 2009 with the poloidal component be- 
coming more dominant. At all epochs the radial magnetic field is predominantly non- 
axisymmetric while the azimuthal field is predominantly axisymmetric with a ring of 
positive azimuthal field around the pole similar to that seen on other active stars. 

Key words: line : profiles - Stars : activity - imaging - magnetic fields - Stars : 
individual : HD 141943 - starspots 



1 INTRODUCTION 

The generation of magnetic fields is arguably one of the 
most important process operating inside a star afi'ecting ev- 
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erything from the angular momentum evolution of the star 
through to the habitability of any planets around the star. 
In the solar case the magnetic dynamo is believed to op- 
erate in an interface layer between the differentially rotat- 
ing convective z one and the radiat ive zone which rotates 
as a solid-body (|Parker et al.l [l993 ). However, for young, 
rapidly-rotating solar-type stars evidence is growing that 
such stars may in fact have a differen t dynamo mechanism. 
It has been suggested (|Donati et al.|[20 03:) that such stars 
may in fact house distributed dynamos, i.e. dynamos which 
operate across the entire convective zone, rather than be- 
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ing restricted to the interface-layer as in the solar case. 
The strongest evidence for this are large regions of near- 
surface a zimuthal magnetic field seen on such stars (i.e. 
iDonati et al... 2003. : , Marsden et al,.2006a l. Such regions are 
beheved to be the near-surface toroidal components of the 
large-scale dynamo field. In a solar-like interface-layer dy- 
namo such fields should not be seen near the stellar surface 
and thus it is postulated that a distributed dynamo is oper- 
ating in young, rapidly-rotating solar-type stars. 

Most current dynamo models are based on our under- 
standing of th e Sun and are tailored to reprodu ce solar ob- 
servations (see lParker et al]| 19931 : [Charbonneau 2005). Such 
models usually involve an interface layer dynamo which, as 
discussed, may not apply to young solar-type stars, but there 
are some different dynamo models that have been exam- 
ined. For example, the dynamo operating in fully convec- 
tive stars (i.e. Browning 2008), or near-surface dynamos (i.e. 
lBranderiburgll2005l) or rapidly-rotatin g solar models with no 
interface layer (i.e. iBrown et al.ll2010l ). Still the operation of 
stellar magnetic dynamos is not well understood. 

One of the most direct ways of observing the stellar 
dynamo is through the observation of the global magnetic 
field on the surface of a star. Zeeman Doppler imaging (ZDI, 
ISemell 1 19891 : iDonati fc Brownl liggTl ) has been used for a 



urations of solar-type stars (i.e. 


Donati et al. 


'l99i 119991: 


Donati & Collier CameronI 1997: 


Donati ,1999 


: Petit et all 


2004albl. 20081: Marsden et al.l 2006a: Jeffers & Donatil 20081: 


Dunstone et ahlboosl) and there is now a 


growing body of 



observations of the surface topologies of other stars that any 
stell ar dynamo models ne e ds to be able to reproduce. 

IDonati fc Landstreetl (|2009l ) have summarised many of 
the observations of the magnetic fields of non-degenerate 
stars and for young solar-type stars (i.e. those with Rossby 
numbers < 1 and more massive than ~0.5 M0) they find 
that these stars produce substantial toroidal fields and have 
mostl y non-axisy mmetric poloidal fields. For more mature 



stars, 



Petit et al. I (I2OO8) have shown that for solar-type stars 
with rotation periods faster than ~12 days the global mag- 
netic field of the star appears to be dominated by toroidal 
field with poloidal field being by far the most dominant field 
configuration for stars with slower rotation rates. This could 
indicate a change in the dynamo mechanism for solar-type 
stars rotating faster than ~12 days. Such observations have 
yet to be modelled in stellar dynamo theory. 

Doppler images of the surface spot topology of young 
solar-type stars show that t hey almost universally ha ve 
large polar spot f e atures (i.e. Baxnes et al. 2000l . l2001al lbl: 



large polar spot t e atures I i.e. [Jjajnes et al.l I200UI . 1200 
IDonati et alll2003l : iMarsden et al.ll2005bl . 120063). This 
contrast to the predi c tions of current dyn a mo theory (i.e . 
ISchussler etld] Il996l : ICranzer et al] I2OO0I : ICranzej 120041 ') 

and it has been suggested that a strong meridional flow may 
be responsible for transporting these spots from their emer- 
gence latitudes to the polar regions. This may also explain 
the mixed polarity of m agnetic fields seen o n the polar re- 
gions of many stars (see lMackav et a"l]|2004l '). 

Although great progress has been made, our knowledge 
of stellar magnetic dynamos is still in its infancy. How the 
operation of the stellar magnetic dynamo depends on basic 
stellar parameters such as age, mass and rotation rate is still 
unknown. 

One of the most obvious effects of the solar dy- 



namo is the reversal of the Sun's global magnetic topol- 
ogy every ~11 years, but so far only two other solar- 
type stars have shown evi dence of glo bal magnetic polar- 
ity reversals, HD 190771 jPetit et al.l i2009'l and Tau Boo 
IDonati et al.ll2008l : iFares et al.ll2009l ). No strong evidence 
for a global polarity reversal has yet been seen on a young 
solar-type star, even though some, for example AB Dor 
iDonati fc Collier Cameroi3il997l : IDonati et alTl 19991 . |2003|) 



have been observed for a number of years. Thus it is still 
unknown if young solar-type stars undergo regular magnetic 
cycles like the Sun or have chaotic cycl es, as indicated by th e 
Calcium HK emission of young stars (|Baliunas et al.iri995l ). 

For young solar-type stars the available spectropolari- 
metric observations are mainly of K-stars (i.e. IDonati et al.l 
l2003h . At present there is only one young early-G star 
for which spectropolarimetric observations have bee n pub- 
lished in a referc cd jo urnal , HD 171488 (Marsd en et al.l 
l2006al : Ijeffers fc Donatill2008l : Ijeffers et al]l2O10i), with pre- 
hminary results for the late- F star HR 1817 (jMengelll2006l : 
iMarsden et al.|[2"006bl . l2010al ) being published as conference 
proceedings. 

As part of a study into the magnetic topology and cy- 
cl es of young late-F/early -G stars this p aper along w i th tha t 
of IMarsden et all (|2010bl . Paper H) and lWaite etal] (|2010h . 
presents spectropolarimetric observations of two early-G 
pre-main sequence (PMS) stars. This paper and Paper H 
deals with the o bservations of the young early-G star HD 
141943, while the lWaite et"aD (|2010l ) paper deals with obser- 
vations of the similarly aged but mor e massive HD 1065 06. 

HD 141943 was first identified bv lWaite etHI l|2005l ') as 
a potential target for spectropol arimetric ob servations. It is 
active (Log(Lx) ~ 30.7 erg s~^. |Cutispot~et al.. 2002) and 
bright (V = 7.9, ICutispoto et al.|l2002l). It also has a very 
strong Lithium line (Ali = 3.3jCutisDoto et al."200d) in- 
dicating its youth (~15 Myrs, ICutispoto et al. 2003) and 
it is what we class as a moderately rapi d rotator with 
a vsim ~38 km s~^ jlCutispoto et al J I2OO3I ) and a period 
of 2.20 ± 0.03 days JCutispoto et al.lll999l ). According to 
iHillenbrand et al.l 'l|2008l ') HD 141943 is also possibly host to 
a ~85 K debri disk. 

This paper (Paper I) describes the evolution in both 
the brightness (reconstructed at 4 epochs) and magnetic 
(reconstructed at 3 epochs) topologies of HD 141943 taken 
at the Anglo- Australian telescope (AAT). A further paper 
on HD 141943 (Paper II) discusses the differential rotation. 
Ha emission and coronal magnetic field maps reconstructed 
from these observations. 



2 OBSERVATIONS AND DATA REDUCTION 

HD141943 was observed at 4 epochs on the 3.9-m AAT, in 
May 2006, March/ April 2007, April 2009 and March/ April 
2010. The May 2006 observations were standard spectro- 
scopic observations taken over a 6 night period using the 
University College London Echelle Spectrograph (UCLES). 
The 2007, 2009 and 2010 observations were spectropolari- 
metric, observing both left- and right-hand circularly po- 
larised light usi ng the SEMPOL (or SEMELPOL) spec- 
tropolarimeter (jSemel. Donati fc Reed Il993l : IDonati et al.l 
200t) visitor instrument in conjunction with UCLES. The 
observations in 2007, 2009 and 2010 were taken over 11, 7 
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and 12 nights respectively. A log of the observations is given 
in Table [T] The rotational phase of the observations was 
determined from the following ephemeris: 



HJD = 2454195.154 + 2.1820, 



(1) 



where HJD is the Heliocentric Julian date of the observation 
and (f) is the rotational phase. 

The SEMPOL spectropolarimeter involves a fibre feed 
from the Cassegrain focus of the AAT to UCLES, where 
the two polarisation states (in this case left- and right-hand 
circular polarisation) are outputted in two fibres to UCLES 
with both polarisation states being recorded simultaneously 
on the detector. The basic construction of the SEMPOL po- 
larimeter consists of a quarter- wave plate (or quarter- wave 
Fresnel rhomb, depending on instrument setup) and a half- 
wave Fresnel rhomb. Spectropolarimetric observations in cir- 
cularly polarised light (Stokes V) consist of a sequence of 
four exposures. Between the exposures the half-wave Fres- 
nel rhomb is rotated between -1-45° and -45°. Thus the po- 
larisation in each output fibre is alternated resulting in the 
removal of spurious polarisation signals from the telescope 
and polarimeter (at least to a first order approximation). 
Due to throughput loses in the spectropolarimeter setup the 
signal-to-noise (S/N) of spectropolarimetric data is usually 
lower than that obtained for straight spectroscopic observa- 
tions. Further information on the op eration of th e SEM POL 
spectropolarimeter can b e found in ISemel et all (|l993l ) and 
iDonati et all (Il997l . l2003l) . 

As each Stokes V observation consists of 4 exposures 
this means that we generally have 4 times the number of 
Stokes I (intensity) observations as we do Stokes V. Note 
however that due to poor weather during some runs there are 
a few sequences on HD 141943 where we were only able to 
obtain a sequence of two exposures (see Table[T|. The Stokes 
V profile can still be extracted from just two exposures, but 
with typically lower S/N. 

The same detector and wavelength setup was used for 
all 4 observing runs. The detector was the EEV2 CCD with 
2048 X 4096 13.5 ^m square pixels. As the EEV2 is larger 
than the unvignetted field of UCLES a smaller window for- 
mat (2048 X 2746 pixels) was used to reduce readout time. 
Using the 31 gr/mm grating, 46 orders (#129 to #84) were 
observed, giving a full wavelength coverage from ~4380 A 
to ~6810 A. For the May 2006 observations the chip was 
binned by 2 in the spectral direction and a 1 arcsec slit was 
used providing a resolution of ~50,000 (i.e. ^6.0 km s~^). 
For the spectropolarimetric observations (in 2007, 2009 and 
2010) there was no binning and the fibre feed provided a 
resolution of ~70,000 (i.e. ~4.3 km s"^). 

All raw frames were reduced into wavelength calibrated 
spectra using the ESpRIT (Echelle Spectra Re duction: an 
Intera ctive Tool) optimal extraction routines of Do nati et al.l 
l|l997l ). As the Zeeman signatures in atomic lines are ex- 
tremely small (typical relative amplitudes of 0.1 per cent 
or less) we have a pplied the techniqu e of Least-Squares De- 
convolution (LSD. lDonati et al.llT997[ ) to the over 2600 pho- 
tospheric spectral lines in each echelle spectrum in order to 
create a single high S/N profile for each observation. The line 
mask used for the LSD was a G2 line list created from the 
Kurucz atomi c database and ATLAS9 atmospheric models 
l|Kurucj|l99l ). 



The peak S/N of the initial Stokes I observations ranged 
from 30 to 260 (depending on observing conditions) while 
the peak S /N for the Stokes V observations (combining 2 or 
4 exposures) was 50 to 330. LSD has been applied to both 
the Stokes I and Stokes V data resulting in S/N values of 
540 to 930 for the Stokes I profiles and 1450 to 11600 for the 
Stokes V profiles. This corresponds to a multiplex gain of ~4 
- 10 for the Stokes I observations and ~35 for the Stokes V 
observations. The multiplex gain for the Stokes I profiles is 
significantly less than that for the Sto kes V profiles because, 
as pointed out bv lDonati et al.l (Il997l ). the technique of LSD 
appears to be not as suited to Stokes I as it is to Stokes V. In 
addition, the calculated S/N of the Stokes I LSD profiles are 
usually somewhat underestimated (see Section | 4.1 II . Furthe r 
info rmation on LSD can be found in [Ponati et al.l (Il997l ) 



and iKochukhov. Ma kaganiu k fc Piskunov (|2010l ). 

In order to correct for wavelength shifts of instrumental 
origin, each spectrum was shifted to match the Stokes I LSD 
profile of the t elluric lines containe d in the spectrum, as has 
been done by iDonati et all l|2003l ) and other authors. This 
reduces the relative radial velocity shifts of the LSD profiles 
to ± ~0.1 km s"^ 



3 FUNDAMENTAL PARAMETERS OF HD 
141943 

In order to determine the surface topologies of HD 141943 
the ZDI code needs to know several basic stellar parameters, 
namely theiisini, radial velocity, inclination of the rotation 
axis to the observer, and the temperature of the photosphere 
and spots on the stellar surface. The?; sin j and radial veloc- 
ity are determined using the ^'^-minimisation technique of 
iBarnes et all (|200G| 1. This is simply using those parameters 
that give the best fit to the dataset (lowest reduced val- 
ues) . 

According to iHillenbrand et all (|2008l ) HD 141943 is at 
a distance of 67 pc, has an effective temperature of 5805 
K and a luminosity of 2.7 Lq. Now these latter two values 
will be slightly higher for an unspotted s tar, but we do not 
know the spottedness of the star when the IHillenbrand et al.l 
(|2008l ) determinations were made. Given the level of spotted- 
ness shown on HD 141943 (see Section Pl.tfl we have assumed 
a spot coverage of 0% to 5% with the most likely coverage 
being 3%. This assumption, along with a spot-photosphere 
te mperature d i fferen ce of ~1900 K (calculated from Fig. 7 
of lBerdvuiz:inall2005l ). gives a photospheric temperature of 
~5850 K, a spot temperature of ~3950 K and an unspotted 
luminosity of -^2.8 Lq. 

As there are no e rror measurements given for the 
[Hillenbrand et al] l|2008l ) measurements we have assumed an 
error of ± 100 K in the photospheric temperature and ± 0.1 
Lq in the luminosity. Thes e, along with the uncerta in spot 
coverage at the time of the IHillenbrand et al] (|2008l ) deter- 
minations and a t;sini of 35 ± 0.5 km s~^, imply that HD 
141943 has a radius of ~1.6 ± 0.15 Rq and an inclination 
angle of ~70°]|^^qo . With a high inclination angle our imag- 
ing code (see Section [4]) has difficulty determining if features 
are located in the Northern or Southern hemisphere so we 
have limited the inclination angle range to 70° ± 10°. Ta- 
ble[2]lists the basic stellar parameters of HD 141943 that we 
have determined for this study. As mentioned, there are no 
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Table 1. Logs of spectroscopic and spectropolarimetric AAT observations of HD 141943 for May 2006, March/April 2007, April 2009 
and March/April 2010. The first two columns give the UT dates and times of the mid point of each observation, while the third column 
gives the exposure time. The fourth column gives the rotational phase of the observations calculated from equation [T] For the May 2006 
observations the rotational phases have had 153.0 added to them, the March/April 2007 rotational phases have had 3.0 added to them, 
the April 2009 rotational phases have had 336.0 subtracted from them, while for the March/April 2010 observations the rotational phases 
have had 497.0 subtracted from them. This is so that all phases are shown as positive. The fifth column gives the signal-to-noise ratio of 
the resultant Stoves V LSD profile (only for the 2007, 2009 and 2010 spectropolarimetric data). 



UT date 


UT time 


Exp. 


Rot. phase 






(sec.) 


{+ 153.0) 


06 May 06 


18:44:03 to 


400, 300, 


0.448 to 




19:46:17 


5x600 


0.468 


07 May 06 


17:48:15 to 


9x600 


0.889 to 




19:17:00 




0.917 


08 May 06 


17:16:24 to 


13x600 


1.337 to 




19:27:28 




1.379 


09 May 06 


17:26:23 to 


13x600 


1.798 to 




19:38:43 




1.840 


10 May 06 


16:53:01 to 


15x600 


2.246 to 




19:30:00 




2.296 


11 May 06 


17:39:12 to 


11x600 


2.719 to 




19:27:05 




2.753 



UT date UT time Exp. Rot. phase S/Nlsd 

(sec.) (- 336.0) 

07 Apr 09 13:31:28 4x600 0.349 3827 

16:16:43 4x600 0.405 3327 

18:45:33 4x600 0.449 2251 

08 Apr 09 13:25:37 4x600 0.805 5216 

18:51:41 4x800 0.909 7225 

09 Apr 09 13:59:27 4x800 1.274 5050 

17:07:59 4x800 1.334 5760 

10 Apr 09 14:33:60 2x800 1.744 5902 
13 Apr 09 13:38:03 2x800 3.101 4730 



UT date 


UT time 


Exp. 

(sec.) 


Rot. phase 
{+ 3.0) 


S/Nlsd 


30 Mar 07 


12:30:43 


4x600 


0.649 


4948 


31 Mar 07 


12:20:25 


4x600 


1.104 


4507 


01 Apr 07 


12:26:09 


4x600 


1.564 


4656 




18:40:53 


4x600 


1.684 


4544 


02 Apr 07 


12:23:33 


4x600 


2.022 


5959 




18:14:40 


4x600 


2.134 


7715 


03 Apr 07 


12:26:18 


4x600 


2.481 


4606 




18:56:19 


4x600 


2.605 


5354 


04 Apr 07 


13:07:40 


4x600 


2.953 


5009 




18:58:10 


4x600 


3.064 


6980 


05 Apr 07 


12:09:39 


4x600 


3.392 


4873 




15:42:02 


4x600 


3.460 


4599 




18:40:21 


4x600 


3.517 


6617 


06 Apr 07 


12:33:12 


4x600 


3.858 


4915 




16:18:35 


4x600 


3.930 


6553 




18:35:54 


4x600 


3.974 


6700 


07 Apr 07 


12:18:45 


4x600 


4.312 


7035 


08 Apr 07 


12:29:38 


4x600 


4.774 


6963 




16:11:27 


4x600 


4.844 


6734 




18:47:16 


4x600 


4.894 


6847 


09 Apr 07 


12:18:25 


4x600 


5.228 


5992 




16:11:12 


4x600 


5.303 


7662 




18:40:21 


4x600 


5.350 


8100 



UT date 


UT time 


Exp. 
(sec.) 


Rot. phase 
(- 497.0) 


S/Nlsd 


25 Mar 


10 


13:50:14 


4x800 


0.674 


10144 






18:57:49 


4x800 


0.772 


6921 


26 Mar 


10 


14:08:19 


4x800 


1.138 


9113 






17:03:17 


4x800 


1.194 


2034 






19:17:19 


4x800 


1.237 


4127 


27 Mar 


10 


14:57:05 


4x800 


1.612 


8385 






16:49:52 


2x800 


1.648 


4045 






17:19:28 


2x800 


1.657 


1470 






19:06:03 


4x800 


1.691 


8869 


28 Mar 


10 


14:33:02 


4x800 


2.063 


8330 






17:05:00 


4x800 


2.111 


5814 


31 Mar 


10 


13:06:47 


4x800 


3.410 


5950 






19:01:29 


4x800 


3.523 


9673 


01 Apr 


10 


13:30:07 


4x800 


3.876 


7274 






18:45:07 


4x800 


3.976 


11575 


02 Apr 


10 


13:03:07 


4x800 


4.326 


9440 






18:45:48 


4x800 


4.435 


8337 


03 Apr 


10 


13:16:00 


4x800 


4.788 


8470 






15:50:33 


4x800 


4.838 


9398 






18:55:30 


4x800 


4.896 


4473 


04 Apr 


10 


14:22:06 


4x800 


5.268 


7227 






18:39:53 


4x800 


5.350 


5186 


05 Apr 


10 


12:05:53 


4x800 


5.683 


4779 






13:05:15 


4x800 


5.702 


5559 



error measurements given for the iHillenbrand et all (j2008l ) 
measurements of effective temperature and luminosity, so 
the estimations of the errors in the stellar parameters of HD 
141943 may well be underestimates. 



From the values for photospheric temperature and 
unspotted luminosity we can determine the mass and 
age of HD 141943 from th e evolutionary models of 
ISiess. Dufour fc Forestiiiil l|2000l ). see Fig.[T] This makes HD 
141943 a ~1.3 M© star of age ^17 Myrs, in reason able agree- 
ment with the values from Cutispoto et al.l |2003h , and thus 
it can be classified as a PMS star. 



3.1 Is HD 141943 a binary? 



iNordstrom et al.l (|2004h took 7 radial velocity measurements 
of HD 141943 over a 2979 day period (epochs and measure- 
ments were not given). Based on the fact that the standard 
error of each measurement was less than the standard devi- 
ation of the measurements they determined that HD 141943 
only has a 31.4 per cent chance of having a constant radial 
velocity. 

Our own radial velocity measurements also show some 
level of varia t ion. T aking the radial ve locity measurement of 
IWaite et all ()2005[ ) along with ours (jCutispoto et al ] |l999l 
has some radial velocity measurements but the errors are 
too large to be useful) we have listed them in Table [3] to 
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3.80 3.75 3.70 3.65 3.60 3.55 

Log(T,ff) 

Figure 1. The evolution ary status of HD 141943 from the pre- 
main sequence models of ISiess et al ] l|200(t ). using the 2007 up- 
dates. The solid lines are pre-main sequence evolutionary tracks 
for stellar masses of 1.0 Mq to 1.5 in 0.1 steps, with the 
bold solid line denoting the 1.0 Mq track. The dashed lines are 
age isochrones for stellar ages of 10 Myr, 15 Myr, 20 Myr, and 
the Zero- Age Main-Sequence (bold dashed line). 



Table 3. Log of radial velocity measurements of HD 141943 from 
these observations plus that of Waito et al. (2005). For the 2006, 
2007, 2009 and 2010 datasets the HJD given is for the mid-point 
of the observations. The errors in the radial velocity measures 
represent a 3ct change in the reduced values (see text). 



Data set 


HJD 


Radial Velocity 




(2450000.0-I-) 


(km s-i) 


2005 


2396.828 


-1.5 ± 1.0 


2006 


3864.801 


-0.6 ± 0.1 


2007 


4195.154 


0.1 ± 0.1 


2009 


4932.071 


0.1 ± 0.1 


2010 


5286.564 


0.2 ± 0.1 



the noise level. Therefore a regularisation scheme is usually 
introduced to determine a unique solution with the role of 
the regularisation scheme being r educed with high quality 
data. Our imaging code uses the ISkilhng fc BrvanI (|l984h 
maximum-entropy regularisation scheme. This scheme has 
the effect of producing images with the minimum amount of 
information required to fit the data to the noise level. 



Table 2. Fundamental parameters of HD 141943 used/found in 
th is study. Th e age and mass are from the theoretical isochrones 
of ISiess et al 

] feooo) given in Fig. [T] and the equatorial rotation 
period has been determined from the differential rotation, see 
Paper H 



Parameter 


value 


Age 


~17 Myrs 


Mass 


~1.3 Mq 


Photospheric temperature 


5850 ± 100 K" 


Spot temperature 


~3950 K 


Unspotted luminosity 


2.8 ± 0.1 Lg, 


Stellar radius 


1.6 ± 0.15 R© 


vsini 


35.0 ± 0.5 km s"! 


Radial velocity (tirad) 


see Table |3] 


Inclination angle (i) 


70° ± 10° 


Equatorial rotation period (Poq) 


~2.182 days 


"assumed error. 



highlight the level of variation. The radial velocity for the 
2006, 2007, 2009 and 2010 datasets has been calculated as 
part of the Doppler imaging process (see Section [4]) and 
are for the entire dataset as a whole. Given that the radial 
velocity for 2007, 2009 and 2010 is constant (within error 
bars) it would mean that if HD 141943 is a binary then it 
would appear to have an elliptical orbit. 



4 IMAGES OF HD 141943 

A surface image of a rapidly-rotating star can be gener- 
ated through the inversion of a time series of Stokes I 
(brightness images) or Stokes V (magnetic images) pro- 
files. The brightness and magnetic topologies of HD 14 1943 
were created using the ZDI code of lBrown et al.l (|l99ll ) and 
iDonati fc Brownl l(l997h using a linear limb-darkening coef- 
ficient of 0.66. As the inversion is an ill-posed problem, an 
infinite number of solutions can be found that fit the data to 



4.1 Brightness images 

The imaging code to reconstruct the spot features on HD 
141943 uses a two-temperature model (one for the spots 
an d the other for th e unsp otted photosphere) as described 
by ICollier CameronI (|l992l ). Each pixel on the stellar sur- 
face is reconstructed for spot occupancy (the local rela- 
tive area occupied by cool spots) and varies from (no 
spots) to 1 (maximum spotted ness) . Following on from 
lUnruh fc Collier CameronI (|l995l ). who showed that there is 
little change in spot topology when using a synthetic pro- 
file over profiles created from slowly rotating comparison 
stars, we have used Gaussian profiles to represent the pro- 
files of both the spots and the photosphere. This is almost 
standard practice now with a number of previous works do- 
ing the same (i.e. iP etit et al.| l2004bl : iMarsden et al. I l2005bl . 
l2006a : .Jeffers fc Donati 2008i ) . As has done for the two pre- 
vious young early-G stars for which spot topologies have 
been obtained by this method (R58, Marsden et al.|[2005bl 
and HD 171488, iMarsden etal] l2006ah we have used the 
same Gaussian to represent both the spot and photosphere. 
The full-width at half-maximum (FWHM) of the Gaussians 
used were slightly different for the spectroscopic and spec- 
tropolarimetric data sets, with a FWHM of 10 km used 
for the spectroscopic data set and a FWHM of 9 km s~^ 
used for the spectropolarimetric data sets. These FWHM 
were determined from matching the FWHM of the Moon's 
LSD profile taken with the same instrumental setup during 
each observing run. 

The maximum-entropy brightness images for HD 
141943 for the 4 observing epochs are shown in Fig. [5] (2006 
observations), the top- left image of Fig. [3] (2007 observa- 
tions), the top-left image of Fig. |4] (2009 observations) and 
the top-left image of Fig [S] (2010 observations). These im- 
ages were created fitting the data down to reduced values 
of 0.2, 0.3, 0.4 and 0.3 for the 2006, 2007, 2009 and 20 10 
datasets respectively. As explained in lPetit et al.l (|2004bl ) a 
reduced value smaller than unity can be achieved be- 
cause, as mentioned in Section (2] the S/N calculated for the 
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Figure 2. Maximum entropy briglitness image reconstruetion for 
HD 141943, May 2006. The image is a flattened polar projection 
looking down on the north pole and extending down to -30° lati- 
tude. The bold line denotes the equator and the dashed lines are 
-|-30° and -1-60° latitude parallels. The radial ticks outside the 
plot indicate the phases at which the star was observed. The im- 
age has a spot filling factor of 0.021 (or 2.1 per cent). The image 
incorporates the measured surface differential rotation, see Paper 
II. 



Stokes I LSD profiles are underestimated. This has no ef- 
fect on the maps produced. Fits of the modelled profiles to 
the observed LSD profiles are given in Fig. [6] (2006 observa- 
tions), Fig. [3 (2007 observations). Fig. [8] (2009 observations) 
and Fig [9] (2010 observations). The epochs of the observa- 
tions, as calculated from the midpoint of each dataset, are: 
2006.352, 2007.257, 2009.273 and 2010.244 in decimal years. 

The brightness images of HD 141943 show that at all 4 
epochs HD 141943 had a smallish polar spot and numerous 
low-latitude features situated mostly between the equator 
and ~-|-30° latitude. In the May 2006 image (Fig. [2| there 
appears to be some features below the equator and above 
-1-30°. However, as there are no observations around these 
phases it could well be that the code has had trouble de- 
termining the latitude of spot features that are only seen in 
the wings of the profiles. 

The spot filling factor (the level of spot coverage over 
the entire stellar surface) is also very similar over the 4 
epochs. In 2006 the spot coverage was 2.1 per cent, in 2007 it 
was 3.1 per cent, in 2009 2.7 per cent and in 2010 it was 2.9 
per cent. The reason for the slightly lower value of spot cov- 
erage for the 2006 dataset may well be due to the limited 
phase coverage of the observations, see Fig. [2l Thus some 
spots in these observing gaps may not have been recovered. 

The variation in spot occupancy with stellar latitude is 
given in Fig. IIOI which plots fractional spottedness versus 
stellar latitude. Fractional spottedness is defined as: 



F{e) 



s{e)cos{e)d6 



(2) 



where, S{9) is the average spot occupancy at latitude and 
dd is the latitude width of each latitude ring. 

It has been reported that active solar-type stars show 
evidence of active longitudes, where certain longitudes (usu- 
ally two longitudes around 180° apart) are more active than 
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Figure 8. Maximum entropy fits for the Stokes I LSD profiles of 
HD 141943, April 2009. As described in Fig.\E\ the dots represent 
the observed LSD profiles while the lines represent the fits to the 
profiles produced by the imaging code. 
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Figure 10. Fractional spottedness versus stellar latitude for 
HD141943. Fractional spottedness is based on the average spot 
occupancy at each latitude and is defined by equation [2] 
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Figure 3. Maximum entropy brightness and magnetic image reconstructions for HD 141943, March/April 2007. The images are flattened 
polar projections as described in Fig. |2] The scale in the magnetic images is in Gauss. The brightness image (top-left) has a spot filling 
factor of 0.031 (or 3.1 per cent), while the magnetic images have a mean field modulus of 91.3 G. The images have been created with 
the inclusion of the surface differential rotation of the star, see Paper II. 



others on the stellar surfa ce (i.e. iBerdvugina fc Tuomineril 
1 19981 : Ijarvinen et al.l [20051 ). In order to test this we have 
plotted the average spottedness versus stellar rotational 
phase for all 4 epochs of HD 141943 observations. This is 
displayed in Fig. 1111 In case the polar spot is aflecting the 
results we have plotted the average spottedness for both (a) 
0° to +90° and (b) 0° to +60° in latitude. 

For the 2006 and 2009 datasets there is very little ev- 
idence for active longitudes on the stellar surface and the 
limited phases observed in these datasets may well be influ- 
encing this result. For the well sampled 2007 dataset, and 
looking at the top-left image in Fig. |3l there appears to 
be a slight enhancement of the lower-latitude spots around 
phases ^0.25 and ~0.55. However as shown in Fig. Illf b) 
these enhancements are rather broad. For the 2010 dataset 
there appears to be an increase in the spot coverage at phase 
~0.00. However, this is reduced when only considering low- 
latitude features (Fig. Illf b)! and when looking at the top- 
left image in Fig.[5]the enhancement could be caused by the 



extension of the polar spot to slightly lower-latitudes around 
phase ~0.00 in the 2010 epoch. There would not appear to 
be any evidence of another active longitude on the star in 
2010. We do not feel that these represent strong evidence 
for active longitudes on HD 141943. 



4.2 Magnetic images 

The magnetic field topology for HD 141943 has been 
rec onstructed for three epochs using the modelling 
of iDonati fc Brownl (|l997l ) and including the spherical 
harmonic ex pansion of the surface magnetic field by 
iDonati et al. 

I |2006). The magnetic reconstruction is done 
assuming a generalised potential field plus a toroidal field 
using a high order {I < 30) spherical harmonic expansion. A 
limit of /max = 30 was chosen for the spherical harmonic ex- 
pansion, as beyond this limit the magnetic topologies remain 
essentially unchanged. 

Various weighting schemes can be applied to the spher- 
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Figure 4. Maximum entropy briglitness and magnetic image reconstructions for HD 141943, April 2009. The images are flattened polar 
projections as described in Fig. |2] The brightness image (top-left) has a spot filling factor of 0.027 (or 2.7 per cent), while the magnetic 
images have a mean field modulus of 36.6 G. The images incorporate the measured surface differential rotation, see Paper II. Note the 
different scale in the magnetic images compared to Fig. |3] 



leal harmonic expansion so that the reconstruction favours 
different magnetic field topologies. Following the principle 
of Occam's razor, we have used a weighting scheme that 
favours "simpler" magnetic fields (i.e. those with lower I 
values) while still reconstructing a similar overall magnetic 
topology to an unweighted reconstruction. Details of the 
spheri cal harmonic technique can be found in lDonati et al.l 
ll2006ll . 

It should be noted that the April 2009 dataset can be 
modelled using a purely potential field without the need to 
include a toroidal field, but this may be due to the limited 
data in the 2009 observations. The March/ April 2007 and 
March/April 2010 datasets however, cannot be modelled to 
a reduced x'^ value of 1.0 without the inclusion of a toroidal 
field component. Thus we have assumed a potential field 
plus toroidal field when modelling all three datasets. 

Using the Stokes V LSD profiles the magnetic imaging 
code reconstructs images of radial, azimuthal and merid- 
ional field on the stellar surface, assuming a weak magnetic 



field and a constant Gauss ian profile over the stellar sur- 
face (see lDonati et all 120031 ). Radial field is defined as field 
in/out of the stellar surface, with positive field being field 
lines directed outward from the surface. Azimuthal field is 
defined as field wrapped around the rotational axis of the 
star, with positive being defined as counterclockwise. Merid- 
ional field is defined as field following lines of longitude north 
and south, with positive being northward. 

The reconstructed magnetic topologies from the 
March/April 2007, April 2009 and March/April 2010 epochs 
are given in Fig. O Fig. [3] and Fig. [S] respectively. The fits 
to the observed Stokes V LSD profiles are given in Fig. [12] 
(2007 observations), Fig.[l3](2009 observations) and Fig. [Til 
(2010 observations). All three epochs were fitted to a re- 
duced value of 1.0. Table [4] list the magnetic quantities 
derived from the magnetic maps in Fig. [3] Fig. [4|and Fig. [5| 
and F ig. [T51 plots these values out. As done in TPetit et alj 
(|2008l ) we have varied the input parameters (namely inclina- 
tion angle, vsmi, rotational period and differential rotation) 
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Figure 5. Maximum entropy brightness and magnetic image reconstructions for HD 141943, March/April 2010. The images are flattened 
polar projections as described in Fig. [2] The brightness image (top-left) has a spot filling factor of 0.029 (or 2.9 per cent), while the 
magnetic images have a mean field modulus of 70.9 G. The images incorporate the measured surface differential rotation, see Paper II. 
Note the different scale in the magnetic images compared to Fig. |3] and Fig|4] 



by the errors given in Table [2] and Paper II to determine 
the possible errors in the magnetic quantities given in Ta- 
ble U These errors apply only for the simplified weighting 
scheme we have used in our reconstructions (see above) and 
assuming a potential plus toroidal field in the magnetic field 
reconstruction. They are thus likely to be lower-limits. 

The reconstructed magnetic topology maps shown in 
Fig. [3l Fig. |4] and Fig. [5] show the maximum entropy recon- 
struction of the large-scale magnetic topology of HD 141943 
at the three epochs, 2007.257, 2009.273 and 2010.244. The 
magnetic imaging code mainly reconstructs the large-scale 
magnetic topology on the stellar surface of the star, as mag- 
netic flux that is contained in small dipoles below the res- 
olution limit of the observations (~11° in longitude at the 
stellar equator) is likely to cancel out and not be recovered. 

As has been previously described (i.e. iDonati fc Brownl 

1 19971 ) the technique of ZDI sometimes suffers cross-talk be- 
tween radial and meridional components of the magnetic 
field for low-latitude features. However, this effect is mainly 



present in stars with low inclination angles. Given the rela- 
tively large inclination angle of HD 141943 {i ~ 70°) this ef- 
fect should be minimal. Conversely the sensitivity of ZDI to 
low-latitude meridional fields decreases significantly with an 
increase in stellar inclination, and thus the meridional field 
topologies of HD 141943 (lower-right images of Fig. O Fig. [J] 
and Fig. [5| may well miss some meridional field that has not 
been reconstructed. However, given the lack of meridional 
field seen on most active stars we believe that if there is any 
missing field it should be minimal. 

In the April 2009 observations, there is a large gap in 
the observations between a phase of ~0.45 to ~0.75. This is 
most likely to have an effect on the recovered radial magnetic 
field (as ZDI is mostly sensitive to radial field around the 
phase of the observation and to azimuthal field around 0.2 
phase from the observation). This could explain the slightly 
lower intensity of the radial magnetic features (top-right im- 
age of Fig. |4]) at these phases and may well mean that the 
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Figure 6. Maximum entropy fits for the Stokes I LSD profiles of HD 141943, May 2006. The dots represent the observed LSD profiles 
while the lines represent the fits to the profiles produced by the imaging code. Each profile is shifted down by 0.01 for graphical purposes. 
The rotational phases at which the observations took place are indicated to the right of each profile. 



Table 4. Magnetic quantities derived from the magnetic maps in Fig.|3l Fig.|4]and Fig. (5] The first column is the epoch of the observation, 
while the second column is the mean field modulus over the stellar surface. The third and fourth columns are the per cent of the large 
scale magnetic energy in the reconstructed poloidal and toroidal field components. The fifth, sixth and seventh columns are the per cent 
of the poloidal magnetic energy in dipole {1 = 1), quadrupole {I = 2) and octupole components (I = Z), while the eighth, ninth and tenth 
columns are the same for the toroidal field. Columns 11 and 12 are the per cent of the poloidal and toroidal field that are axisymmotric 
(m = 0). 
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per cent of poloidal magnetic energy in Table |4] is slightly 
underestimated for the April 2009 observations. 

The first thing to notice from Table |4] is that the mean 
magnetic field of the 2007 and 2010 maps is significantly 
higher than that of the 2009 maps (~91 G and ~71 G com- 
pared to ~37 G). This could either be a significant change 
in the magnetic field strength on HD 141943 during 2009 
or it may well be due to the limited amount (and poorer 
quality) of observations in the 2009 dataset compared with 



that of 2007 and 2010. To test this we took only a limited 
number of profiles from the 2007 and 2010 datasets that had 
a similar phase coverage to that of the 2009 dataset. This 
had the result of decreasing the mean magnetic field to ~66 
G and ~51 G for the 2007 and 2010 datasets respectively. In 
addition, we artificially increased the S/N of the Apr 2009 
profiles by fitting the data to a reduced of 0.8, which had 
the effect of almost doubling the mean magnetic field on 
the star to ~60 G. Thus we believe that the poor quality of 
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Figure 7. Maximum entropy fits for tfie Stolces I LSD profiles of HD 141943, Marcli/ April 2007. Again tlie dots represent tlie observed 
LSD profiles while the lines represent the fits to the profiles produced by the imaging code as described in Fig. (6] 



the April 2009 data could well be the reason for the change 
in mean magnetic field on the star, rather than any actual 
change in the magnetic field strength of the star. 

The other obvious feature of Table U is the observation 
that the balance between poloidal and toroidal magnetic 
energy changed dramatically in 2009 compared to the 2007 
and 2010 epochs, being roughly balanced in 2007 and 2010 
but being dominated by the poloidal component in 2009. 
Using the magnetic maps created from a limited number of 
profiles in 2007 and 2010 (see above) we tested to see if the 
limited dataset in April 2009 could have affected the ratio 
of poloidal and toroidal magnetic energy seen in the 2009 
observations. Reducing the number of profiles used in the 
creation of the magnetic map had very little effect on the 
balance between poloidal and toroidal field in both the 2007 
and 2010 datasets (still almost evenly balanced). Thus we 
believe that the change in the ratio of poloidal and toroidal 
magnetic energy in the 2009 epoch may well be real, but we 
cannot rule out the poor dataset in 2009 being responsible. 

The major effect that limiting the number of observa- 
tions in the 2007 and 2010 datasets had was to slightly de- 
crease the complexity of the magnetic field (with more en- 
ergy in octupole components or lower) and this may explain 



the slight lowering of the complexity of the poloidal field 
seen in the 2009 dataset. 

As was done for the spot features, the variation in mag- 
netic field with stellar latitude is given in Fig. [TS] This is 
again determined using equation [2j but with S{9) represent- 
ing the value of the magnetic field at each latitude rather 
than spot occupancy. We have also determined the varia- 
tion in magnetic field with longitude, averaging the absolute 
value of the magnetic field at each longitude. This is shown 
in Fig. [T71 

Fig. [TH] shows that the radial magnetic field of HD 
141943 is predominantly negative at all epochs for latitudes 
above ~0° . The azimuthal field is predominantly positive at 
all latitudes in 2007 and 2010 but in 2009 it is positive at 
high latitudes and predominantly negative at latitudes be- 
low ^--^60° . Although we feel that this change in the azimuthal 
field, along with the change in the ratio of poloidal/toroidal 
field, is possible evidence for a changing magnetic field on 
HD 141943, we cannot rule out the poor 2009 dataset being 
responsible. As with the brightness images. Fig. [T7] shows 
little or no evidence of active longitudes in the magnetic 
fields on HD 141943. 
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Figure 9. Maximum entropy fits for the Stokes I LSD profiles of HD 141943, March/April 2010. As described in Fig.|6l tiie dots represent 
tlie observed LSD profiles while the lines represent the fits to the profiles produced by the imaging code. 
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Figure 15. Plot of the magnetic parameters of HD 141943 given in Table [4] for the poloidal field (right-hand plot) and toroidal field 
(left-hand plot). The solid line is the per cent of the large scale energy in the reconstructed poloidal or toroidal field components, the 
percent of the poloidal or toroidal magnetic energy in dipole (dot-dashed line), quadrupole (dashed line) and octupole (dotted line) 
components are also given. The percent of the poloidal or toroidal field that is axisymmetric is plotted as a triple-dot-dashed line. 
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Figure 17. Average of the absolute value of the radial magnetic field (upper plots) and the azimuthal magnetic field (lower plots) versus 
stellar rotational phase for HD141943, averaged over 0° to +90° (left-hand plots) and 0° to +60° (right-hand plots) latitude. As was 
done in Fig. 1111 the phase has been binned into 0.1 steps to remove any small scale variations. 



5 DISCUSSION 

HD 141943 is only the se cond (or thir d including the re- 
sults for HD 106506 by Wai te et al.koiol ') young early-G star 
for which the large-scale magnetic to pology has been deter- 
mined, the other beirig HD 171488 (iMarsden et all l2006al : 
IJeffers fc Donatil l2008l : Ijeffers et all I2OI0I) . However, there 
have been five young early-G stars for which spot maps and 
differential rotation measures have been determined (and a 
number of others that have just spot maps, see Table 4 in 
IStrassmeier et al.ll2003l ). Along with HD 141 943, HD 171488 
and H D 106506, the ot her two stars are, R 58 (|Marsden et al.l 
l2005a lH) and LQ Lup (jDonati et al.ll2000l ). For comparative 
purposes in this discussion the stellar parameters of all these 
five stars are given in Table [S] 



5.1 Spot topology 

The surface spot topology of HD 141943 at four epochs is 
shown in Fig. [5] Fig. [3] (top-left image) , Fig. U (top-left im- 
age) and Fig. [5] (top-left image) . The 4 maps show that the 
spot topology of HD 141943 was remarkably consistent over 
the span of the observations (~4 years). All maps show that 
HD 141943 has a smallish polar spot with a number of lower- 
latitude features situated predominantly between the equa- 
tor and +30° latitude, with only the 2010 epoch showing 
some significant spot features between +30° and +60° lati- 
tude at around phase ~0.90. The total spot coverage for all 



4 epochs is also very similar ranging from 2.1 per cent to 3.1 
per cent. 

Comparing these maps to those of other young early- 
G st ars created using the same im aging code, such a s 



R58 llMarsden et all l2005b|) LQ Lup (iDonati et al.l 



2000l) . 



l2008l : 



HD 17148 8 ('Marsden et al.' '2006a'; 'jeffers & Donati 
IJeffers et al . 2010) and HD 106506 (Waite et a l. 201(J), it is 
quite noticeable that while all stars appear to have some 
lower-latitude features, the polar spot on HD 141943 is sig- 
nificantly smal l er than that shown by these other targets. 
iMarsden et al. I ('2005b') has shown that starspot mapping as- 
suming high stellar inclination angles leads to a dramatic in- 
crease in the amount of polar spot features needed to match 
the observed deviations in the LSD profiles. Given that HD 
141943 has a higher stellar inclination than the other stars 
listed in Table [5] it is not an incorrectly determined stellar 
inclination that is responsible for the small polar spot on 
HD 141943. 

As can be seen in Table O HD 141943 is the second 
youngest and second most massive (behind HD 106506) of 
the young early-G stars so far imaged using this code. It 
also has the shallowest convective zone (as a function of the 
stellar radius) and is the slowest rotator of the five. Which of 
these factors is responsible for the small polar spot seen on 
HD 141943 is still open for debate. The young early-K star 
LQ Hya has a rotation rate somewhere between that of HD 
141943 and HD 171488 (P ~ 1.6 d) and has a convective zone 
depth of ^^0.29 R*. In some spot maps LQ Hya also show a 
smallish polar spot (jDonati et al.ll2003l ). although the maps 
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Table 5. Comparison of the stellar parameters of the five young early-G stars that have had their spot topology and surface differential 
rotation measured using Dopplcr imaging. The first three, HD 141943, HD 106506 and HD 171488, have also had their magnetic 
topo logies im aged. Except where noted, the data for HD 141943 are taken from this paper and Paper II, th e data for HD 106506 are 
froml Waite~et al . (2010 ) and that for HD 171488 are from lStrassmeier et ahl J2003 ) and M arsden et alj l l2006al '). The data for R58 comes 
from M arsden et al.l lj2005bl ) and the data for LQ Lup from lDonatr'et''anil2000l ) . The values for the depth of the convective zone are 
from lSiess et al.l tiOOff) . 



Parameter 



HD 141943 



HD 106506 



HD 171488 



R58 



LQ Lup 



(B-V) 0.65" 0.605 0.62" 

Age (Myrs) ~17 <10 30 - 50 

Mass (Mq) ~1.3 1.5 ± 0.1 1.20 ± 0.02 

Radius (R©) 1.6 ± 0.1 2.15 ± 0.26 1.15 ± 0.08 

Inclination {°) 70 ± 10 65 ± 5 60 ± 10 

Convective zone (R*) ~0.16 [0.26 Rq] ~0.22 [0.47 Rq] ~0.21 [0.24 R©] 

Equatorial rotation period (d) ~2.2 ~1.4 ~1.3 



o.ei*" 

35 ± 5 
1.15 ± 0.05 

60 ± 10 
-0.21 [0.25 R©] 
~0.56 



0.69'= 
25 ± 10 
1.16 ± 0.04 
1.22 ± 0.12 
35 ± 5 
-0.26 [0.32 Kq] 
-0.31 



"from lCutispoto et al] l l2002tl : ''from lRandichI l l200ll ). dereddened value; '^from lWichmann et ahl l ll997ll . 
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Figure 11. Average spottedness versus stellar rotational phase 
for HD141943, averaged over (a) 0° to +90° and (b) 0° to +60° 
latitude. The phase has been binned into 0.1 steps to remove any 
small scale variations. 



that show this are often those th at do not have well s ampled 
observational phases (i.e. Fig. 13. lDonati et al.ll2003h . In ad- 
dition, prelim inary observations of the late-F star HR 1817 
l|Marsden et al.. .2006b'l. with a rotation period of ~1.0 d and 
a convective zone depth of ~0.17 R*, have shown it to have 
a small, but intense, polar spot. A significantly larger sam- 
ple size is required before any firm conclusions can be made 
about the relationship between polar spot size/intensity and 
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Figure 12. Maximum entropy fits for the Stokes V LSD profiles 
of HD 141943, March/April 2007. The dots represent the observed 
LSD profiles while the lines represent the fits to the profiles pro- 
duced by the imaging code. Each profile is shifted down by 0.003 
for graphical purposes. The rotational phases at which the obser- 
vations took place are indicated to the right of each profile while 
the 1(7 error bars are given to the left of each profile. 
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Figure 13. Maximum entropy fits for tiie Stoltes V LSD profiles 
of HD 141943, April 2009. As explained in Fig.[l2l the dots rep- 
resent the observed LSD profiles while the lines represent the fits 
to the profiles produced by the imaging code. 



basic stellar parameters. Indeed, given the poss ible variation 
in th e size of the polar spot seen on LQ Hya jDonati et al.l 
120031 ) it may well be that young active stars go through 
cycles where the size of th e polar spot varies as has been 
me ntioned by Barnes! (|2005l ). 

iGranzei l|2004l ) has used thin flux tube models to deter- 
mine the latitude distribution of spot emergence for young 
rapidly-rotating stars. Although these models are based on 
a solar-like interface dynamo (which may not be applicable 
to such stars, see Section 15. 2p it is still interesting to see 
how well the models compare to the results for HD 141943. 
Fig. 1 101 shows that at all 4 epochs HD 141943 has a peak in 
fractional spottedness around 0° to 4-40° latitude, with an- 
other peak around+70° to+85°. For a PMS star of 1.0 Mq 
the models of iGranzeil (|2004| ) predict that the star should 
have a spot emergence distribution of between +30° to +70° 
latitude with a peak around +50°. For a 1.7 M© PMS star 
this distribution is more spread out with a range from +20° 
to +70° with a peak around +30° and +40°. The distribu- 
tion of fractional spottedness shown in Fig. [10] is not that 
dissimilar to the comb ination o f the 1.0 Mq and 1.7 Mq 
distributions given by iGranze 

2 (|2004l ) with the exception 
that the higher-latitude peak in spottedness is at a higher 
latitude in our images than predicted by the models. How- 
ever, the models predict that for the more Zero- Age Main- 
Sequence (ZAMS) age stars, such as most of the other stars 
in Table O the latitude of peak distribution in spots is pre- 
dicted to decrease, but the more intense polar spots of these 
stars would appear to indicate that this is not happening. 
As suggested bv lGranzeil (|2004l ) and others, it is possible 
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Figure 14. Maximum entropy fits for the Stokes V LSD profiles 
of HD 141943, March/ April 2010. As explained in Fig. [12] the 
dots represent the observed LSD profiles while the lines represent 
the fits to the profiles produced by the imaging code. 



a strong meridional flow is pushing most of the spots, which 
form at lower-latitudes, to the polar region, with stars with 
deeper convective zones possibly having stronger meridional 
flows and thus forming more intense polar spots. Such a 
mechanism might well explain the bimodal spot distribu- 
tion seen in Fig. [10] with the lower peak in the distribution 
representing the latitude at which the flux tubes emerge 
and the polar spot those spots that have been driven to 
the poles by the meridional flow. Such a mechanism may 
also explain the non-uniform polar spots seen on a num- 
ber of stars, if the polar spot is being formed by a number 
of smaller spot feature s. This is still just speculative, but 
I Weber. Str assmcier & Washuettll (2005) have shown tenta- 
tive evidence for large poleward meridional flows on early-K 
giants. 



5.2 Large-scale magnetic topology 

The maps of the large-scale magnetic topology on HD 
141943 from three epochs (Fig. [3] Fig. [4] and Fig. [5]) show 
that the radial field on HD 141943 had a mixed polarity 
at all latitudes for all epochs. The azimuthal field appears 
to be dominated by a ring of positive field around the pole 
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Figure 16. Fractional magnetic field for (a) the radial magnetic 
field and (b) the azimuthal magnetic field versus stellar latitude 
for HD141943. This is based on the absolute value of the magnetic 
field at each latitude and is defined by equation [2] with S{9) now 
representing the value of the magnetic field. 



at all epochs, although the intensity of the polar ring is re- 
duced in 2009. This is shown in the results from Table [l] 
which show that the toroidal field is predominately axisym- 
metric (~70 - 75 per cent) while the poloidal field is pre- 
dominately non-axisymmetric (~70 - 80 per cent), with the 
ratio not appearing to change significantly between the three 
epochs. The other obvious feature of the azimuthal field is 
the increase in the amount of negative field on the stellar 
surface in the 2009 image compared to the 2007 and 2010 
images. As these surface regions of radial and azimuthal field 
are believed to be the poloidal and toroidal components of 
the large-scale dynamo field, the latitudinal distribution and 
evolution of these regions are an important window onto the 
operation of the stellar magnetic dynamo. 

There appears to be no large-scale magnetic polarity re- 
versal on HD 141943 over the 3 years of obser vations. This is 



simila r to the young early-G st a r HD 171488 ( Marsden et al.l 
l2006al : I JeflFers Donatii l2008l : IJeflFers et all l2010ll and the 



youn g early-K stars AB Dor and LQ Hva (|Donati et al] 
l2003l ) none of which have shown evidence of magnetic polar- 
ity reversals over several years of observations. In contrast, 
a magnetic polarity reve rsal has been see n on the mature 
Sun-hke star HD 190771 (jPetit et aLlbOOgl ) in 3 years of ob- 
servations and the mature late-F star Tau Boo has shown 
two magnetic polarity reversals in the space of ^2 years 



(jDonati et al.ll2008l : [Fares et aLlbOOglV Howev er, Tau Boo is 
host to a "Hot Jupiter" ( Butler et al.l Il997l l that may be 
affecting the star's magnetic cycle. 

An axisymmetric field in the form of a ring of azimuthal 
field around the pole appears to be a common element of 
active solar-type stars with such a ring also seen on all re- 
constructed images of the young early - G stars HD 171488 
(Marsden et al. 2006a'; JefiFers & Donati 2008'; 'jefiFe rs et all 
,2010,) and HD 106506 (,Waite et al. 2010 1. as weU as being 
seen on a number of im ages of the young early-K star AB 
Dor l|Donati et al.ll2003h . Such a ring has also been seen on 
more evolved so lar-type stars such as the early-G FK Com 
star HD 199178 (iPetit et al."2004lj ) and the early-K R SCVn 
star HR 1099 (|Donati et al...2003i : iPetit et al.|[2004al '). Such 
large-scale regions of azimuthal magnetic field near the sur- 
face of these active stars has lead to the belief that such 
stars have a fundamentally different dynamo in operation 
than that of the solar-type interface-layer dynamo, as such 
regions should not be seen near the stellar surface if a solar 
dynamo were in action. It is thought that such stars con- 
tain a distributed dynamo, one that operates throughout 
the stellar convective zone rather than b eing restricted to 
the interface layer (see lDonati et af]|2003l l. 

Such a co ncept has rec e ntly b een supported by theoret- 
ical models of iBrown et~al] (|201G| ') which have shown that a 
stable magnetic dynamo can exist without an interface-layer. 
These models are based on a rapidly-rotating Sun (up to 3 
times the current solar rotation rate) and also show large 
regions of near-surface "wreaths" of longitudinal magnetic 
field. These wreaths are similar to the rings of azimuthal 
field seen on active solar-type stars, except that they are 
at lower latitudes. Perhaps the higher Coriolis force experi- 
enced by the more rapidly rotating stars such as HD 141943 
(rotating around 10 times the solar value) are pushing these 
wreaths to higher latitudes. 

Table 3] shows that in all epochs both the poloidal and 
toroidal magnetic fields are very complex with over 50 per 
cent of the magnetic energy in orders higher than an oc- 
tupole (except for the poloidal magnetic energy in 2009 
which has only 45 per cent of the magnetic energy in orders 
higher than an octupole). There also appears to be evolution 
in the magnetic field between the epochs. In 2007 the mag- 
netic energy is reasonably evenly balanced between poloidal 
and toroidal components (47 per cent to 52 per cent respec- 
tively) while in 2009 the poloidal field is extremely dominant 
(82 per cent to 17 per cent) and in 2010 the magnetic energy 
is again balanced between the poloidal and toroidal compo- 
nents (50 per cent in each). As explained in Section [4.21 the 
dataset in 2009 is significantly poorer than those in 2007 
and 2010 and this may be having an influence on the recon- 
struction of the magnetic topology for this epoch. However, 
when we tried to reconstruct the magnetic topology for the 
2007 and 2010 epochs using a small subset of the data simi- 
lar to that of the 2009 dataset there was virtually no change 
in the ratio of poloidal and toroidal field. Thus we feel that 
the change in the poloidal/toroidal ratio seen in 2009 may 
well be real, but unfortunately we cannot prove this with 
any certainty. 

iDonati fc Landstreetl (|2009| ) has looked at the changing 
magnetic field topologies on stars of different masses a nd ro- 
tation rates (see Fig. 3 in lDonati fc Landstreetl |2009| ). This 
shows that stars around the mass and rotation rate of HD 
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141943 (actually preliminary results from the March/ April 
2007 dataset on HD 141943 are plotted in the figure) fall into 
an 'active' dynamo regime, which includes stars with Rossby 
number < 1 and more massive than ~0.5 Mp,. These stars 
produce predominantly non-axisymmetric poloidal magnetic 
field topologies as well as having significant amounts of 
toroidal field. At all three epochs studied in this paper 
HD 141943 does appear to have a predominantly non- 
axisymmetric poloidal field (and a predominantly axisym- 
metric toroidal field), but, as mentioned above, the ratio of 
poloidal to toroidal magnetic field appears to change for the 
2009 epoch. 

A change in the ratio of poloidal to toroidal magnetic 
field has been seen on another star. The previously men- 
tioned mature late-F star Tau Boo (period ~3 days) shows 
large changes in its ratio of poloidal/toroidal field during 
a magnetic cycle, with toroidal field dominating before a 

J lobal polarity sw itch and poloidal field dominating after 
Fares et al.ll2009l '). This could indicate that the changes seen 
in the ratio of poloidal and toroidal field on HD 141943 are 
part of a stellar magnetic cycle on the star, however it would 
mean that it has undergone an extremely rapid evolution of 
its magnetic field, but has yet to undergo a magnetic po- 
larity reversal. It would be useful to follow HD 141943 on 
at least yearly timescales to see whether the magnetic field 
does in fact undergo rapid magnetic evolution and polarity 
reversals, such as that evidenced by Tau Boo. 

As part of the Sun's magnetic cycle, the ratio of dipole 
to quadrupole components of the magnetic field changes. As 
can been seen in Table |4] there are significant changes in the 
percent of magnetic energy in dipole, quadrupole and oc- 
tupole components for both the poloidal and toroidal fields. 
For example the per cent of poloidal magnetic energy in the 
dipole component changes from ~6 per cent in 2007 to ~29 
per cent in 2009 and down to ~6 per cent in 2010. The most 
significant changes appear to occur for the 2009 dataset. 
Assuming that the change seen in the magnetic topology of 
HD 141943 in 2009 is real, a comparison of the total mag- 
netic energy (both poloidal and toroidal components) of HD 
141943 appears to show that HD 141943 has also undergone 
a change in the complexity of its magnetic field from 2007 
to 2009 and then back again in 2010. From Tableglthe total 
dipole component of the magnetic energy stays reasonably 
constant over all three epochs (~24 per cent in 2007, ~25 
per cent in 2009 and ~22 per cent in 2010). The same can be 
said for the total octupole component (~9 per cent in 2007, 
~7 per cent in 2009 and ~12 per cent in 2010). However, the 
total quadrupole component appears to change markedly in 
2009 compared to the other two epochs (~7 per cent in 2007, 
~18 per cent in 2009 and ~9 per cent in 2010) mostly due 
to an increase in the quadrupole component of the poloidal 
field. Now, as explained in Section 14.21 the poorer observa- 
tions in 2009 may play a role in this change but this could 
be a change imparted by magnetic evolution on HD 141943. 



6 CONCLUSIONS 

HD 141943 is one of the youngest and most massive stars 
for which we have analysed the surface magnetic topologies 
for. For young stars it has the shallowest convective zone 
that we have yet studied. We have presented reconstructed 



brightness (at 4 epochs) and magnetic (at 3 epochs) topolo- 
gies of HD 141943. During the four year timebase of the 
brightness images the spot distribution on HD 141943 has 
changed very little with a smallish polar spot and a num- 
ber of lower-latitude features at around 0° to -f 30° latitude, 
seen at all epochs. 

Over three years of observations the pattern of the mag- 
netic topology of HD 141943 also looks at first glance to 
have experienced relatively little change, with positive and 
negative radial magnetic field seen at all latitudes and a 
ring of positive azimuthal magnetic field seen at all epochs. 
At all epochs the large-scale poloidal magnetic field of HD 
141943 is mostly non-axisymmetric while the toroidal field 
is predominantly axisymmetric. The reconstructed magnetic 
topologies are rather complex with over 50 per cent of the 
magnetic energy in components higher than an octupole. 
When the magnetic images were analysed in more detail 
we found tentative evidence for a change in the magnetic 
field in 2009. The ratio of poloidal to toroidal field on HD 
141943 goes from almost balanced in 2007 to being heav- 
ily dominated by poloidal magnetic field in 2009 and back 
to balanced in 2010. If real, this variation would indicate 
magnetic evolution on HD 141943. 
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